Josephson junctions as threshold detectors for full counting statistics 
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We discuss how threshold detectors can be used for a direct measurement of the full distribution 
of current fluctuations and how to exploit Josephson junctions in this respect. We propose a scheme 
to characterize the full counting statistics (FCS) from the current dependence of the escape rate 
measured. We illustrate the scheme with explicit results for tunnel, diffusive and quasi-ballistic 
mesoscopic conductors. 



Quantum noise in electron transport is an actively de- 
veloping field. Noise measurements provide exclusive in- 
formation about microscopic mechanisms of the trans- 
port that can hardly be obtained by other means 0, • 
Still, the experiments in the field neither match the in- 
tensive theoretical development nor gather all informa- 
tion about electric fluctuation. Indeed, the concept of 
full counting statistics pioneered in Q allows one to pre- 
dict the non-Gaussian distribution function of the current 
measured during a time interval r, P T (I). This distri- 
bution is characterized by an infinite set of cumulants 
<C I n !S>. A traditional noise measurement only assesses 
the second cumulant of this set discarding the rest. Re- 
cent pioneering work reports a successful measurement 
of the third cumulant but there is a long way to go if 
one measured the cumulants one by one. It would be ad- 
vantageous to measure the distribution function directly 
and thus to get all cumulants at once, thereby collecting 
the wealth of information being currently discarded. 

Why is such a measurement difficult? The probabili- 
ties to measure correspond to big deviations of the cur- 
rent from its average value, |7 — (7}| ~ (7), and are there- 
fore exponentially small. For instance, in the shot noise 
regime P T (7) ~ exp {-(I)G{I / (I))r/e), 0(1/ (I)) ~ 1 
being the function to characterize. One has to concen- 
trate on very rare measurement outcomes that occur with 
probability exp(— (7)r/e) w 0. Such measurements can 
only be carried out with threshold detectors that discrim- 
inate these rare events. Let us discuss an ideal thresh- 
old detector that measures the current during the time 
interval r, and gives a signal if the current measured ex- 
ceeds the threshold current Ith- The signal probability 
will then be proportional to P T (7 t h). To give a realistic 
illustration, a detector that measures a tunnel junction 
with (7) = 10 pA in the time interval r = 10~ 6 s would 
go off once an hour if 7 t h = 2(7) and once in 10~ 4 s if 
7 t h = 1-5(7). Therefore, if one measures the rate of the 
detector signals as a function of 7 t h, one directly assesses 
the full counting statistics. 

Albeit realistic detectors are not ideal. There are three 
important factors that can either hinder the interpreta- 
tion of such a measurement or even prevent the measure- 
ment: (i) a realistic detector hardly measures the current 
averaged over a certain time interval r. It is dispersive, 



being usually more sensitive to longer and smaller cur- 
rent fluctuations rather than to bigger and shorter ones, 
(ii) The detector may produce a significant feedback on 
the system measured when it starts to signal, thereby dis- 
rupting its noise properties, (iii) The detector could just 
go off by itself, for instance, due to quantum tunneling. 

A Josephson junction seems to be a natural threshold 
detector for current fluctuations. It can be viewed as a 
particle in a washboard potential , the superconducting 
phase difference <f> across the junction corresponding to 
the particle's coordinate. The junction is in zero-voltage 
state provided the current does not exceed the critical 
value corresponding to the critical tilt of the washboard 
potential. <p is trapped in one of the minima of the poten- 
tial, which is separated by a barrier from the neighbor- 
ing one. A current fluctuation that exceeds the critical 
threshold sets </> into motion and the junction gives a 
signal — a voltage pulse that lasts till <j) is retrapped in a 
different minimum. 

In this paper we address the feasibility of Josephson 
junction systems for measuring the full distribution of 
current fluctuations in a mesoscopic conductor. Our re- 
sults are as follows. The Josephson junction is a realistic 
detector, all three factors mentioned are in play. Albeit 
one can measure the distribution provided the width of 
the barrier <fio ^> 1 . This can be realized by a flux division 
using two inductances. Under these conditions, the third 
factor is of no importance and the first and second factor 
do not hinder the unambiguous correspondence between 
FCS and the escape rate of the junction as a function 
of 7 t h- These theoretical results open the way to direct 
experimental observation of FCS. 

The circuit under consideration consists of a normal 
coherent conductor with conductance G in series with 
the Josephson junction(system) (Fig. 1). The system 
is biased with voltage source V 3> k^T/e. This assures 
that the normal conductor is in the shot noise regime 
and its actual temperature is not relevant. In addition, 
we inject extra current 7b that controls the slope of the 
Josephson washboard potential. 

If fluctuations are neglected, this system can be de- 
scribed with the celebrated model of resistively shunted 
junction 0. The normal conductor is a source of 
non-Gaussian current fluctuations that instantly tilt the 
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FIG. 1: A voltage biased mesoscopic conductor with conduc- 
tance G provides the noise source for a threshold detector 
which is characterized by its threshold current J t h. Ih is an 
additional current bias. The inset shows a possible realiza- 
tion of the detector with two inductances £1,2, a Josephson 
junction and a negative inductance element. 



washboard potential and can lead to an escape of <fi from 
the minimum. The escape gives rise to an observable 
voltage pulse. The escape rate in the same or similar 
systems has been studied for a variety of noise sources 
and potentials 0, 0> B @ • To our knowlegde, the non- 
Gaussian noise sources that are characterized by FCS 
were not adressed yet. 

To proceed, we begin with the fully quantum mechani- 
cal description of the system in terms of a Keldysh action 
for a single variable <p ll(l | , that incorporates informa- 
tion about FCS of the normal conductor and the proper- 
ties of the Josephson junction. We calculate the escape 
rate by considering saddle-point trajectories of the ac- 
tion, A, that connect the potential minimum with the 
nearest potential maximum. With exponential accuracy, 
the rate is given by T ~ exp(— ImA/h). 

The action consists of two terms, A = Aj + A-^, cor- 
responding to the elements of the circuit. We denote 
by (jr 1 the phases on the forward/backward parts of the 
Keldysh contour and also use symmetrized combinations 
of these 4>,x = (<£i = ' = < /'~)/2- The junction part reads in 
a standard way 



Aj = dtl U(q>+(t)) 



h 2 C 



^ +2 (t) -{0 



C being the self-capacitance of the junction, U(cj>) be- 
ing the Josephson energy with the current bias term 
included: —U{<j>) — (h/2e)(I c cos + Ib<t>) for a single 
junction. Further we concentrate on overdamped junc- 
tions where C -C G 2 Tij '(2eJ c ) and neglect the capacitance 
term. The normal conductor part we write following [ll) 
in quasi-stationary approximation 

AN = Ye G J *(^-^(*))S(x(*)); (2) 

where S characterizes the FCS and the preceding fac- 
tor is just the voltage drop over the normal conductor. 
The distribution of current noise is given by the Fourier 



transform of S. Derivatives with respect to x generate 
the moments of the distribution. 

A coherent conductor can be presented by a set of 
transmission eigenvalues T n and S is given by Levitov's 
formula yl 



^(x) = ^E ln ( 1+T »( elx - 1 )) 



G 



(3) 



Gq being the conductance quantum. Concrete forms of 
S(x) for specific conductors will be given below. At x —> 
0, S can be expanded in x, S ~ ix ~ X 2 F/2> F being the 
Fano factor that describes the suppression of shot noise 
in comparison with the Poisson value [l| . 

This quasi-stationary approximation is only valid if 
the typical time r of the motion along the saddle-point 
trajectory is long in comparison with h/eV, that is, 
eVr 3> h. To check the validity of this, we precede the 
results with simple qualitative estimations. 

Let us consider an arbitrary barrier with the width 
0o and height Uq ~ (7i/e)/ t h0o- The detection time 
can be estimated equating the potential energy term 
and the term with 0, G<f>o(h/e) x/ T — Uqx/<Po, X be- 
ing a typical value along the trajectory. This gives 
t ~ (h/ eV)(f>o(If / Ith)- The quasi-stationary approxi- 
mation thus holds provided If = GV » Ith/<f>o- Let 
us estimate x by equating the term which is quadratic 
in x an d the potential term. This gives x — hh/If if 
Ith *C if, X — 1 otherwise. We see that if 0o < 1 then 
^ « 1. The latter implies that S(x) can be expanded 
near x — and only the first two cumulants are relevant: 
no chance to see the effect of FCS. However, if <fio 1, 
X can become of the order of unity without violating the 
quasi-stationary approximation, and one can observe the 
FCS. The quasi-stationary approximation remains valid 
for x < 0o- 

The resulting rate can be estimated as logT ~ 
4> {G/Gq)x- If 0o lj this reduces to logT ~ 
4>Q{G/Go)It\ l /If. In the opposite limit, the estimation 
for the rate reads logT ~ cfio (G/G Q)H(/ t h /It), 2 being 
a dimensionless function ~ 1. It is important to note 
that these expressions match the quantum tunneling rate 



r±\ logV h ~ U T/h ~ {G/Gq)4>1 provided eVr ~ h. There- 



fore the quasi-stationary approximation is valid when the 
quantum tunneling rate is negligible and the third factor 
mentioned in the introduction is not relevant. For equi- 
librium systems, the situation corresponds to the well- 
known crossover between thermally activated and quan- 
tum processes at k^Tr ~ h ||. 

We proceed with the quantitative solution. The trajec- 
tories we are looking at start at t — > —00 in the minimum 
of the potential with <f> = m i n ,X = an d approach the 
maximum = max ,X = at t — > 00. They obey the 
equations of motion 



ox 



u(4, + (t)) u(<t>-(t)) + Y e G{v - ^(t))s( x (t))(h 



3 



o = |[W + W)-f(« f ))] + i(i) 3 «i|. 

It is important to note that these equations have a simple 
integral of motion 



i{u{<f> + ) - u {<f>-)) + -hS(x)=i 



(6) 



1 = for saddle-point trajectories of interest. The full 
action along the trajectory then reads 



2e 2 

k 2 g 



A 



dt<fiS( X ) 



d^sum (7) 



where in the last relation \ is expressed in terms of <f> by 
means of EqUJ 

Let us start with the results for 4>q ~ 1. In this case, 
one expands the action in terms of \ keeping terms of 
the first and second order only. This immediately yields 
X = i4e(dU / d<fi) / (hF If) . The general answer for the 
escape rate can be obtained at any shape of the barrier 
and reads: 



exp 



knT* 



k B T* = eVF/2 (8) 



This is thermal activation with an effective temperature 
given by the noise in the normal conductor. A similar 
effect of noise was envisaged in a recent article f° r 
the phase diffusion regime. 

How to realize a device where the barrier width <pQ 3> 
1? It can not be just a single Josephson junction since 
the phase drop on the junction can not exceed w. We can 
make the phase drop over the junction much smaller than 
the phase drop over the device by flux division with two 
inductances in series provided h/(el c ) < L± <C L2 
(see inset Fig. 1). However, this is not enough since the 
energy of the device would be dominated by that of the 
inductances, cx 4> 2 /(L2 + Li). This parabolic background 
shall be compensated with a negative inductance — {L\ + 
L2) in parallel. Such negative inductance can be made 
with the aid of an active circ uit [l3l.ll4| or properly biased 
Josephson junction system 15] . This provides a wide 
barrier U{4>). 

We notice that any function U(<j>) can be approximated 
by a cubic parabola if the tilting of the washboard poten- 
tial is close to the critical value. This is why we choose 
the cubic parabola form 



dU 



2e 



hh 



(9) 



for actual calculations. It is convenient to require that 
the barrier does not change if we change If. This can be 
done by a corresponding change of 1^. To simplify this 
further, we notice that x ^ 4>o so that 



dU 



^ubstitution into Eq. HJ gives <fi in terms of x 




hh V iX 



(11) 



Combining this with Eq|7| we obtain the escape rates as 
a function of Jth/Jf f° r anv given FCS. 

To stress similarities and differences with thermal ac- 
tivation, we present the results in the form of Arrhenius- 
like plots. We plot logT in units of (G/Gq)4>o ver- 
sus the dimensionless I t h/I{- Thermal activation with 
the effective temperature given by ||SJ| would give a 
straight line (dashed lines in the plot). By virtue of 
our approach, the rates should exceed the quantum limit 
logT^ ~ (G/Gq)<Pq. This means that the rates should 
saturate at this value provided If — ► 0. For each choice 
of S(x) we pl°t two curves corresponding to two possible 
signs of V with respect to the current via the junction. 
For forward bias, the barrier is crossed when the fluctu- 
ating current is smaller than the average current. For 
backward bias, the barrier is crossed if the fluctuating 
current is bigger than the average value. The difference 
between two curves thus reflects the asymmetry of the 
current distribution with respect to the average current. 

In Fig. 2, left panel, we present the results for a tun- 
nel junction (St(x) — e%x ~~ 1) an d a diffusive conductor 
(S d (x) = (l/4)arccosh 2 (2e ix - 1)) 0. All curves ap- 
proach the dashed thermal activation lines at If ^S> Ith- 
Since the tunnel junction is more noisy (F = 1 versus 
F = 1/3 for a diffusive conductor), it generally provides 
higher escape rates. However, the difference in functional 
form of the rates remains pronounced even upon rescal- 
ing with factor 3. The most pronounced feature of the 
backward bias curves is a plateau at If — > Ith with subse- 
quent drop to very small escape rates ~ Tn (beyond the 
vertical scale of the plot). This is because the current 
distribution is restricted: shot noise current is always of 
the same sign as the average current. 

A quasi-ballistic conductor presents two peculiari- 
ties of this kind. We choose the transmissions of 




(10) 
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FIG. 2: Escape rates versus Ith/Ii f° r a tunnel (t), diffusive 
(d) and ballistic (b) mesoscopic conductor. refers to 

forward/backward bias respectively. Dashed lines correspond 
to the rates due to Gaussian noise. 
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FIG. 3: Optimal current fluctuations (left panel) and feedback 
voltage (right panel) versus time for different conductors and 
different values of I t ^/If. Each line corresponds to one point 
on the curves in figure 2. Note the different voltage scales. 
The labels stand for (branch|7 t h//f): (a) b+|l/5, (b) b+|l/10, 
(c) t-|l/3, (d) d-|l/3, (e) d-|2/3, (f) t-|2/3. 



all channels to be the same, To = 0.8, Sb(x) = 
(1/T ) In (l + T (e 4X - !))• In this case, the current dis- 
tribution is restricted from both sides: the maximum cur- 
rent fluctuation can not exceed the ballistic limit I\ = 
If /To. From this we conclude that the barrier can not 
be crossed at forward bias if 7th > (1/Ib — l)If = 0.25/f, 
as seen in the right panel of Fig. 2. The rate becomes 
increasingly smaller upon approaching this threshold. 

There is an unambiguous correspondence between the 
rates as a function of if/ith and S, that is given by Eqs. 
0, Ijllfl and can be used to characterize the FCS from 
the rates measured. However, this relation is implicit and 
more complicated than that of an ideal detector. Appar- 
ently, this complication is due to the first and second 
factor mentioned in the introduction. To look at it in 
more detail, we compute the optimal current and voltage 
fluctuations that switch the detector. 

The optimal current fluctuations are plotted in the 
left panel of Fig. 3 for different conductors and If. The 
curves are symmetric owing to the symmetry of the cubic 
parabola potential. Common features are that they all 
reach the threshold current at maximum and their time 
spread is of the order of r. Still, the spread, shape, and 
most importantly, the integral of the current over time, 
varies significantly from curve to curve. This proves that 
the detector in use is dispersive and suffers from the first 
factor mentioned in the introduction. 

The third factor is also in play. When <f> moves, cross- 
ing the potential barrier, the resulting voltage changes 
the voltage drop over the normal conductor thereby af- 
fecting the current fluctuations in there. This feedback 
voltage Vfb is negative for forward bias and positive for 
negative one. We see from the evolution equations that 



Vft_ 
V 



H _ s( x (t)) d x 



2eV 



X(t) dS(x(t)) 



- 1, 



(12) 



so the change in the voltage drop across the junction is 
quite significant if % ~ 1. We check that the negative 
feedback can never change the sign of the voltage for 
S(x) in use. The right panel of Fig. 3 presents voltage 



fluctuations corresponding to the current fluctuations on 
the left panel. Interestingly, the positive feedback can 
be very big on the plateau at the backward bias (curves 
c, f). In this case, the detector seeks to optimize the 
rare fluctuation where almost no current is flowing in the 
normal conductor. The probability of such fluctuations 
is increased upon increasing the voltage drop over the 
conductor so that the detector provides the extra volt- 
age required. Eventually, the feedback can be reduced 
with an extra resistive shunt over the Josephson junction. 
However, this would decrease r and reduce the region of 
applicability of our results. 

To conclude, we proved that Josephson junctions can 
be used as threshold detectors for non-Gaussian noise 
produced by coherent conductors. Our theoretical re- 
sults facilitate a new type of electric noise measurement: 
direct measurement of the full distribution of transferred 
charge. 
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